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We study a minimal flavor violating extension of the MSSM, where higher dimensional operators 
in the Kahler potential induce tree level non-holomorphic Higgs couplings that are controlled by 
the scale of the physics beyond the MSSM and analyze their possible impact on CP violation in 
B s and Bd mixing. We consider results on the time dependent CP asymmetries in B s — > ipcf> and 
B s — > ipfo from LHCb, in B s — > iptfi from CDF and DO and in B — > tpKs from the B factories as 
well as the measurement of an anomalous like-sign dimuon charge asymmetry at DO. Taking into 
account the stringent bounds on the branching ratio of the rare B s — ¥ decay, we investigate 

to which extent the framework allows to address the observed (2 — 3)<r discrepancies in fits of the 
unitarity triangle. We find that a non-standard Bd mixing phase, that is in agreement with the 
current bounds on CP violation in B 3 mixing, requires the presence of higher dimensional operators 
both in the Kahler potential and the superpotential. The corresponding region of parameter space 
is characterized by small tan/3 ~ 5, a light Higgs spectrum with masses below < 400 GeV and will 
be probed by future measurements at LHCb. 



INTRODUCTION 

Low Energy observables sensitive to CP Vio- 
lation in b — > s transitions constitute excellent 
probes of possible new CP violating phases in 
extensions of the Standard Model (SM). Indeed, 
as CP Violation in such processes is predicted 
to be tiny in the SM, evidence for sizable CP 
violation in b — > s transitions would be a clear 
hint for the presence of New Physics (NP). Ex- 
amples of such low energy probes are observ- 
ables that are sensitive to the B s mixing phase, 
as the semi-leptonic asymmetry aj in decays of 
B s mesons to "wrong sign leptons" or the time 
dependent CP asymmetries in the B s — > ip<f> and 
B s -> ipfo decays. 

In the context of generic two Higgs doublet 
models with Minimal Flavor Violation (MFV), 
where the CKM matrix is the only source of 
flavor violation, a large B s mixing phase can 
be realized if additional CP violating phases are 
allowed [1-4] . Simultaneously, these models can 
also address tensions in fits of the CKM matrix 
that seem to indicate a sizable NP contribution 



to the Bd mixing phase [5, 6]. Possible relations 
between a non-standard B s mixing phase and 
the baryon asymmetry of the universe in these 
models have been studied in [7, 8]. 

In the context of the Minimal Supersymmet- 
ric Standard Model (MSSM) on the other hand, 
a MFV soft sector is not sufficient to gener- 
ate sizable NP phases in meson mixing [9, 10] 
due to the strong experimental bound on the 
B s — > fj, + fi~ branching ratio. Non SM-like B s 
and Bd mixing phases in the MSSM require 
new sources of flavor violation in addition to 
the CKM matrix (see e.g. [10, 11] for studies 
of such frameworks). Supersymmetric models 
with MFV can generate large B mixing phases if 
they allow for a strongly reduced muon Yukawa 
coupling such that the BR(B S — > /i + /i~) con- 
straint can be avoided. Such a situation can 
be realized for example in the so-called uplifted 
SUSY Higgs region [12], even though this frame- 
work is strongly constrained by other B physics 
observables and (g — 2)„ [13]. As studied in 
detail in [14], non-negligible corrections to CP 
violating observables in meson mixing can also 
be generated if the MSSM with MFV is ex- 
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tended by the two leading higher dimensional 
operators in the Higgs sector with complex co- 
efficients [15]. Still, the stringent bounds on 
BR(_B S — > only allow for a B s mixing 

phase of S^<p < 0.15 in specific regions of pa- 
rameter space of this model. 

In this work we analyze an extension of the 
MSSM, introducing higher dimensional opera- 
tors not exclusively in the Higgs sector, but 
also considering dimension 5 operators that in- 
duce non-holomorphic Higgs-fermion couplings 
already at the tree level [16]. Such operators 
are a possible source of flavor and CP viola- 
tion. Assuming that these operators obey the 
Minimal Flavor Violation ansatz [17-19], we ex- 
plore to which extent the considered framework 
allows for non-standard CP violation in B mix- 
ing without being in conflict with the bounds on 
the BR(£? S — > For a study of a similar 

framework, that however does not consider CP 
violation in B mixing, see [20]. 

CP VIOLATION IN B MESON MIXING 

The B q mixing amplitude 

M« = Mf 2 - l -T\ 2 , (1) 

consists on an dispersive part, Mf 2 , and an 
absorptive part, T\ 2 . The absorptive part is 
dominated by tree level SM contributions and 
therefore hardly affected in many NP models. 
Throughout this work we will assume that Y\ 2 
has no significant NP contributions. The dis- 
persive part M\ 2 on the other hand is highly 
sensitive to new heavy degrees of freedom. The 
effects of NP in M\ 2 can be parametrized by 

Ml 2 = C q e<\Ml 2 ) SM . (2) 

The main impact of the parameters Cd and C s 
is on the mass differences AMd and AM S , re- 
spectively. The NP phases </>^ p and 4^ p affect 
observables that are sensitive to CP violation in 
B meson mixing, like the semileptonic asymme- 
tries ag L and etg L as well as the time-dependent 



CP asymmetries in Bd — > ipKs, B s — > ipifi 
and B s — > ipfo- In order to constrain the NP 
parameters through measurements of these ob- 
servables, knowledge of the respective SM con- 
tributions is required. As several of the ob- 
servables, in particular AMd and S$k s ~ the 
time dependent CP asymmetry in Bd — > ipKs - 
are a key ingredient in the determination of the 
Unitarity Triangle (UT), a simultaneous fit of 
the CKM parameters and the NP parameters 
as performed in [5] is the most consistent ap- 
proach. Here we focus mainly on the impact of 
the recent improvements on the determination 
of the B s mixing phase at CDF and DO [21, 22] 
and in particular at LHCb [23]. As these mea- 
surements have no significant effect on the de- 
termination of the CKM parameters and the NP 
parameters other than <^ p , we consider a sim- 
plified approach and take the CKM parameters 
as well as Cd and C s from the generic NP fit 
in [5], and fit only the NP parameters <^ p and 
fid P - We expect this approach to give a good es- 
timate of the allowed region of parameter space 
compatible with the present experimental data 
on B meson mixing. In particular we will use 

p = Atg[(V tb V t * d )/(V cb V: d )} = (27.2j:&) (3) 
Ps = Arg [(V tb V t * s )/{V cb V: s )\ = (-1.3 ± 0.1)° 

and the following 2a bounds 

0.62 < C d < 1.15 , 0.79 < C a < 1.23 . (4) 

We now give expressions for the observables 
that are sensitive to CP violation in B meson 
mixing. For the semileptonic asymmetries we 
obtain [5, 24] 

10 4 a^ L ~ (55.2 sin ^ p - 4.80cos^ p )/C d (5) 
10 4 a SL ~ (49. 7 sin (^ p + 0.19 cos <^ P )/C S (6) 

where the uncertainties on the numerical coef- 
ficients are at the level of ~ 15%. For the time 
dependent CP asymmetries one has 

S^K 8 = sin(2/3 + <$ p ) , (7) 
S w =Sv, /0 =sin(2|&|-O. (8) 
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These expressions hold under the usual assump- 
tion that the Bd — > ipKs, B s — > tfxp and 
B s — > ipf decays are dominated by the SM tree 
level amplitudes. Using (3) the corresponding 
SM predictions are 

S SM s = 0.82+0.02 ; s sm = 0M6t o.o^ (9) 

The DO collaboration measured the like-sign 
dimuon charge asymmetry that is predicted to 
be composed out of the semi-leptonic asymme- 
tries in the Bd and B s decays [25] 

4 L = 0.59(4 + 0.414 L 

= (-78.7 ± 19.6)10~ 4 . (10) 

The corresponding SM prediction ^4g L (SM) = 
(-2.8t°j)10~ 4 [5] is roughly a factor 25 be- 
low the central value in (10) and differs from 
it by 3.9<7. The value (10) updates an earlier 
DO study [26] that found a 3.2<r evidence for 
an anomalous like-sign dimuon charge asymme- 
try. A separate extraction of the semileptonic 
asymmetries results in [25] 

(4(,D0) = (-12±52)10- 4 , (11) 
a s SL (D0) = (-181 ± 106)10~ 4 . (12) 

The results in Eqs. (10), (11) and (12) hint to- 
wards large negative values for the NP phases 
4>^ p and, in particular, <^ p . 

Interestingly enough, there is a (2-3)<r tension 
between the SM prediction of S^k s (7) and its 
experimental value [27] 

S"ks =0.67 ±0.02 , (13) 

that is largely driven by the B — > rv measure- 
ments that prefer a large value of \V u b\. This 
tension points in the same direction for the NP 
phase as the data on the like sign dimuon 
charge asymmetry. 

A small preference for a negative NP phase 
in B s mixing was also observed in CDF and DO 
data on the time dependent CP asymmetry in 
B s -> ip<t> that give @ 95% C.L. [21, 22] 

- 1.36 <^ p -2|ft|< 0.26 (CDF), (14) 

— 1.65 <4^ F — 2|/3 S |< 0.24 (DO) . (15) 
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FIG. 1. Allowed ranges for the NP phases 4>^ F and 
at the 1 and 2 a level, taking into account the 
measurements of the time-dependent CP asymme- 
tries in B s -> ip(j} at CDF, DO and LHCb [21, 23], in 
B a -> ipf at LHCb [23] and in B d -¥ ipK s at the 
B factories [27]. The measurement of the like-sign 
dimuon charge asymmetry at DO [25] is included in 
the black dotted contours but not in the red solid 
contours. 

Combining the results from [21] with [26], global 
fits to the data found a B s mixing phase <^ p = 
O(-l) [5, 28]. 

Recently however, LHCb presented results on 
the time dependent CP asymmetries in B s — > 
ipcj> and B s — > tpfo that are consistent with the 
tiny SM prediction and that strongly restrict 
the possible values for a NP phase in B s mix- 
ing [23] 

^ p -2|&|= 0.03 ±0.16 ±0.07. (16) 

In Fig. 1 we show the result of a simple fit of 
the NP phases </>^ p and a^ p to the combined 
LHCb result on the time-dependent CP asym- 
metries in B s —> ip(j> and B s — > ipfo [23], the 
results on the time-dependent CP asymmetry 
in B s -> ip(j) from CDF and DO [21, 22] as well 
as the measurement of S^,k s at the B facto- 
ries [27], using the values for j3 and /3 S , Cd and 
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C s from above. The allowed region is mainly 
determined by the measurements of S^, S^f 
at LHCb and S$Ks a t the B factories, while 
the measurements of at CDF and DO lead 
to a small shift of the central value of (j)^ p to- 
wards a small negative value. We stress that 
the very large value of the like-sign dimuon 
charge asymmetry observed by DO cannot be 
explained given the current data on the time 
dependent CP asymmetries in B s — > ip<j) and 
B s — > tpfo- Using our fit results we find a cen- 
tral value of Ag L = — 11 • 10~ 4 and a 2a range 
of -18 < 10 4 A b SL < -2. This differs from the 
measured value (10) by ~ 3a. Including the 
like-sign dimuon charge asymmetry directly into 
the fit leads only to a small shift towards slightly 
larger negative <^ p values. In the following we 
focus on the fit that does not include the 
measurement. 

Due to the small discrepancy between the ex- 
perimental determination of S^,k s an d its SM 
prediction coming from the UT fits [5, 6], the 
NP phase in Bd mixing shows preference to- 
wards a negative value 

= -0.20t°; ° , (17) 

that is roughly 2a below 0. While also for the 
NP phase in B s mixing we find a slight prefer- 
ence for a small negative value, </>^ p is perfectly 
consistent with zero 

^np = -0.10 ±0.15 . (18) 

Presently, this still leaves room for NP contri- 
butions, but the bound (18) will improve signif- 
icantly in the near future with more data from 
LHCb. 

THE B q -)• n+fi- DECAYS 

The CDF collaboration reported a small ex- 
cess in B s —> candidates [29], leading to 

BR(B S -> m + aOcdf = (1.8±J:|) x lO" 8 , (19) 



No excess has been observed by LHCb and CMS 
that report the following bounds [30] 

BR(B S -> M + /OLHCb < 1-5 • 10- 8 , (20) 

br{b s ->m + /Ocms <i-9-nr 8 . (21) 

Combining the bounds from LHCb and CMS 
one finds [31] 

BR(B S -> m + aOlhc < 1-1 x 10- 8 (22) 

that is only a factor of 3.5 above the SM pre- 
diction [32] 

BR(B S ->■ m + aOsm = (3.2±0.2) x 10~ 9 . (23) 

The current bounds on the Bd — > /i + fJT branch- 
ing ratio [30] are still a factor 40-50 above the 
SM expectation and therefore Bd — > is 
much less constraining than B s — > in 
models with MFV. 

Given the strong bound on BR(£> S — » 
possible neutral Higgs contributions to B mix- 
ing in the MSSM with MFV are strongly con- 
strained [33, 34]. Also in the BMSSM model 
considered in [14] non-standard B mixing phases 
are rather restricted (S$$ < 0.15) and can only 
be generated in particular corners of parameter 
space. 

In the following we present an extension of 
the MSSM that respects the MFV principle 
but allows nonetheless for sizable NP phases 
in B mixing without being in conflict with the 
BR(£> S — > n + fj,~) constraint. 



EFFECTIVE HIGGS - FERMION 
COUPLINGS BEYOND THE MSSM 

At the tree level, the MSSM is a 2 Higgs 
doublet model of type II and the couplings of 
the neutral Higgs bosons to fermions are fla- 
vor conserving. At the loop level on the other 
hand, non-holomorphic Higgs couplings are gen- 
erated and have important consequences. Loop 
induced couplings of down-type quarks and 
charged leptons to the up-type Higgs can lead 
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to large threshold corrections to the correspond- 
ing masses [35] and modify significantly CKM 
matrix elements [36] as well as charged Higgs 
couplings to quarks [37] . Finally, they also gen- 
erate flavor changing neutral Higgs couplings 
that can have a profound impact on flavor phe- 
nomenology [33, 34, 38]. All these effects be- 
come relevant for large values of tan that can 
compensate for the 1-loop suppression. 

We now consider possible extensions of the 
MSSM with new degrees of freedom at a scale 
M ~ several TeV. As long as the SUSY break- 
ing scale ms of the new degrees of freedom is 
small compared to M and as long as M is suf- 
ficiently larger than the scale of the MSSM de- 
grees of freedom, one can describe the effects 
of the Beyond the MSSM (BMSSM) physics 
by higher dimensional operators suppressed by 
1/M [15, 39]. An analysis up to order 1/M 
captures the physics of several MSSM UV ex- 
tensions while the effective description of oth- 
ers need to include 1/M 2 effects [40, 41]. In 
this work we restrict ourselves to the 1/M 
level. We consider both the leading higher di- 
mensional superpotential operators that involve 
only Higgs fields [15] 



breaking F-term 



C D 



2M 



J d 2 6(l + aZ)(H u H d 



(24) 



and in particular also 1/M suppressed 
Kahler potential operators that induce non- 
holomorphic Higgs-fermion couplings already 
at the tree level [15, 16] 1 

CD ~ J d 4 6(l + Z + Z^ + ZZ^) (25) 

x {\ u H\QU + \ d HlQD + XtHtLE) + h.c. . 

In the above expressions, Z is an auxiliary 
dimensionless spurion that develops a SUSY 



1 Possible UV completions that lead to the operators 
in (24) and (25) are discussed e.g. in [15, 16, 20, 40] 
and briefly presented in the Appendix. 



Z ->■ m s 6 2 



(26) 



The phenomenological consequences of the op- 
erators in (24) have been thoroughly studied 
in the literature [14, 15, 40, 42, 43]. In par- 
ticular, they can significantly enhance the tree 
level mass of the lightest Higgs boson of the 
MSSM and lead to sizable mass splittings be- 
tween the two heavy neutral Higgs bosons and 
also the charged Higgs boson. Possible phases 
of the coefficients a and w lead also to scalar- 
pseudoscalar mixing. Their impact in the con- 
text of Higgs and flavor phenomenology has 
been analyzed in [14]. 

The Kahler potential operators (25) modify 
the interactions of Higgs bosons, (s)quarks and 
(s)leptons of the MSSM at the 1/M level [15, 
16, 20]. The supersymmetric part of (25) for 
example leads to corrections of the holomorphic 
MSSM Yukawa couplings 



1 M 1 



(27) 



After SUSY breaking also non-holomorphic 
Higgs-quark couplings are generated 

+^<M)ii HlUEj +h.c. , (28) 

where we now made flavor indices explicit. 
These terms can lead to flavor changing neu- 
tral Higgs vertices and correspondingly to tree 
level contributions to FCNC processes like B 
mixing and B s — > In the following we 

will focus on them. The MFV hypothesis as 
formulated in [19] amounts to the assumption 
that the SU(3) 3 quark flavor symmetry of the 
gauge sector is broken by only two spurions y u 
and yd that transform as (3,3,1) and (3,1,3) 
respectively. Correspondingly, the couplings Y q 
and \ q can be expanded in powers of these spu- 
rions. To keep notation simple and concise, we 
conveniently choose y q = Y q — Y q + j^X q to be 
these spurions. Any other linear combination 
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of Y q and X q leads to equivalent results, 
the non-holomorphic Higgs couplings X q to the 
down quarks one then has 



For responding Yukawa couplings y q 



X d = e Q Y d + £l Y d Y^Y d + e 2 Y u Y^Y d + 
+ e 3 Y u Y^Y d Y^Y d + e^Y d YjY u Y^Y d 



(29) 



where for simplicity we dropped the prime on 
the corrected Yukawa couplings. An analo- 
gous expression holds for the up quark coupling 
which is however not relevant for the following 
discussion. The coefficients £j are generically 
of 0(1) and complex. For later convenience we 
define 





m s 




£o = 


~M £ ° ' 






m s 2 


m S 2 2 


£l = 








m s 2 


ms 2 2 


£~2 = 


w y t e2 , 


£4 = —y t y b s4 



£5 — £0 + £1 
£~6 — £0 + £1 



£2 
£4 



- £3 + £4 



(30) 



The couplings in (29) modify the relation be- 
tween the down quark masses m q and the cor- 



VdV _ VsV _ tan/3 

m d m s 1 + £0 tan j3 
ijbv tan /3 

nib 1 + £5 tan j3 ' 



(31) 



with the Higgs vev v — 174 GeV and we only 
show the leading term in a tan j3 expansion. 

Similar to the quark masses, also the CKM 
matrix receives tan /3 enhanced corrections. The 
relations between the affected elements of the 
bare CKM matrix V° and the physical CKM 
matrix V read (i = 1, 2) 



Yi 

Vu 



■10 



1 + £5 tan j3 
1 + eg tan /3 



(32) 



Finally we also give explicit expressions for 
the corrected flavor changing couplings of right 
handed down quarks with the Higgs bosons. 
The leading tan/3 enhanced terms read (i =/= j) 



C D di™^V*V tj Xi S 4 (c a H - s a h + iA) + ul™^- 



Vii Zij d J R H + 



h.c. 



(33) 



t0 



1 + £ 6 t 

(£2 4 



(£2 







£4) ti 



(l + e 5 tp)(l + e 6 tp) 



1 + e 6 t p 1 + eetp e* 2 + e* 3 (ex 



£3)^/3 



1 + e t p 1 + eltp e 2 + £4 



1 + £q£/3 



(34) 
(35) 



The flavor changing bn — > d l L couplings Xib are 
generated by the £2 and £3 terms and they are 
proportional to nib. The flavor changing — > 
d l R couplings Xbi on the other hand are gener- 
ated by the £2 and £4 terms and suppressed by 
light quark masses. The expression (35) gener- 
alizes the results given in [19, 34, 44] and, to the 
best of our knowledge, has not been presented 



in the literature. 

In the MSSM, the 1 factors can only be 
loop induced. Gluino-down squark loops gen- 
erate for example eo, while e 2 is generated by 
chargino-stop loops. Due to the loop suppres- 
sion, the e factors are generically of 0(0.01) 
in the MSSM. Correspondingly, the corrections 
in (31), (32), (34) and (35) become important 
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only for large values of the ratio of the two Higgs 
vevs v u /vd = tan/3 = tp. On the other hand, 
in generic 2 Higgs doublet models with MFV 
as analyzed in [1, 2], where the e are free pa- 
rameters, moderate values of tan /3 ~ 5 — 10 are 
sufficient to generate 0(1) effects. The same 
is true in the supersymmetric framework con- 
sidered here as long as the BMSSM scale that 
controls the size of the e factors is not too high, 
i.e. m s /M ~ 0.1. 



HIGGS EFFECTS IN FLAVOR PHYSICS 

Observables that are highly sensitive to a 
non-standard Higgs sector are observables in 
meson mixing as well as the branching ratios of 
the rare decays B s ^ — > [i + /jT that receive tree 
level contributions from flavor changing neutral 
Higgs exchange. Charged Higgs effects are rele- 
vant in the B — > X s j decay as well as in the 
B — > tv, B — > Dtv and K — > fiv decays. 
The first one is modified only at the loop level 
while the others receive contributions from tree 
level charged Higgs exchange. However these 
tree level decays turn out to give only mild con- 
straints in regions of parameter space with non- 
standard B meson mixing phases and we do not 
discuss them in detail here although they are 
included in our numerical analysis. For a recent 
study of B — > tv in the context of multi Higgs 
doublet models with MFV see [451 . 



B s - B s Mixing 

The flavor changing neutral Higgs couplings 
in (33) give rise to tree level contributions to 
B s mixing mediated by neutral Higgs exchange. 
These contributions can be described by the fol- 
lowing effective Hamiltonian 

n cS = C 2 {b R s L f + C A {l R s L ){b L s R ) , (36) 



Ca — x sh xi s (V tb V t * s ) 



m b m s 



v- 



Ml 



(38) 



Analogous contributions to B c i mixing can be 
obtained through the replacements s — > d. 2 In 
writing (37) and (38), we assume the decoupling 
limit Ma 3> Mz and treat the effect of the oper- 
ators in (24) in a mass insertion approximation. 
In our numerical analysis instead we work with 
Higgs mass eigenstates that we derive from the 
full Higgs potential including MSSM 2-loop cor- 
rections [46]. 

The Wilson coefficient C4 is proportional to 
m b m q . Consequently it can only lead to sizable 
effects in B s mixing while its impact on Bd mix- 
ing is rather restricted. On the other hand, G 2 
is proportional to m 2 and therefore leads to NP 
contributions of the same size and phase both in 
B s and B^ mixing. We stress that C4 is complex 
only if higher orders of the bottom Yukawa are 
considered in the expansion (29) (see also [47]). 
Indeed, one easily checks that switching off E\, 
£3 and £4 leads to a real C4. The Wilson co- 
efficient C*2 however is also highly sensitive to 
the phases of £~o, £2, w and a. From (37) it 
is clear that Ci is only relevant in presence of 
the higher dimensional operators (24) in the su- 
perpotential and for small Higgs masses not far 
above the electroweak scale. 



The B 3 



Decay 



The B s fi + n~ decay receives tree level con- 
tributions from flavor changing neutral Higgs 
exchange. One finds 

BR(B S -> fi + n~) 



Rb s uu — 



BR(B S 



\A\' + \1-A\ 



)SM 

2 



oti) 4Y (x t ) (l + r e tp) Ml 



(39) 



(40) 



with the Wilson coefficients 



Co 



1 ^2 (v T/ *N2 m b atjj m s 



M M\ 



(37) 



2 Neutral Higgs contributions to Kaon mixing are pro- 
portional to m 2 or m B md in the considered framework 
and therefore negligibly small. 
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with the SM loop function given by Yo(xt) — 
0.96. The above expression assumes again the 
decoupling limit Ma ^> Mz- For small Higgs 
masses, corrections at the 1/M level become 
important and are included in our numerical 
analysis. In writing (40) we also assume that 
the non-holomorphic lepton-Higgs coupling Xe 
is proportional to the lepton Yukawa 



I^A* = ^eeYe = e t Y t . (41) 
M M 



As the NP contribution (40) to the B s -t (i + [i~ 
amplitude grows with tan 3 /3, the large tan j3 
regime of the MSSM is strongly constrained by 
the experimental bound (22). However, as al- 
ready stressed above, the non-holomorphic tree 
level Higgs-fermion couplings allow to generate 
NP contributions to B mixing already for mod- 
erate values of tan/3, where the bound from 
BR(£> S — > /i + /j,~) is considerably relaxed as long 
as the muon Yukawa coupling is not largely en- 
hanced by the tan /3 resummation factors. 



NUMERICAL ANALYSIS 

To obtain the B q mixing amplitudes, we use 
2-loop renormalization group running for the 
Wilson coefficients [48] and the hadronic ma- 
trix elements from [49] . We check compatibility 
of the model with various constraints. 

(i) Vacuum stability: The higher dimensional 
operators in the superpotential can lead to a 
second minimum in the Higgs potential. Re- 
quiring that the electroweak minimum is stable, 
gives a lower bound on the charged Higgs mass 
for given values of the fj, term and the SUSY 
breaking scale ms [14, 50]. 

(ii) Electroweak precision observables can 
constrain regions of parameter space where the 
dimension 5 superpotential operators lead to a 
very heavy SM like Higgs boson or a large split- 
ting between the heavy Higgs bosons. We im- 
plement the S and T parameter following [14]. 
Also the Zbb coupling can be modified signifi- 
cantly by Higgs loops [7, 51]. However the Zbb 



constraint can be avoided if the Higgs-top cou- 
pling is suppressed by non-holomorphic correc- 
tions and therefore we do not include it in the 
numerical analysis. 

(iii) Electric Dipole Moments (EDMs) can be 
induced both by phases of the higher dimen- 
sional operators in the Kahler potential and 
the superpotential. Experimentally accessible 
EDMs, like the EDMs of Thallium, Mercury or 
the neutron can be generated by the electron 
and quark EDMs and chromo EDMs (CEDMs), 
by CP violating 4 fermion operators [52] as well 
as by the Weinberg 3 gluon operator [53] . 

The fermion (C)EDMs are generated at the 
1-loop level by sparticle loops that are sen- 
sitive to the phase of the Higgs vev [14, 43] 
as well as to possible complex 1/M correc- 
tions to the sfermion mass matrices that are 
induced by modified Higgs-sfermion couplings 
after electroweak symmetry breaking. The lat- 
ter corrections can however always be avoided 
if the parameters £o an d sg are real. At 
the 2-loop level, Barr-Zee diagrams contribute 
to the (C)EDMs [2, 54, 55]. They are di- 
rectly sensitive to both the phases in the non- 
holomorphic Higgs couplings and the scalar- 
pseudoscalar mixing in the Higgs sector. 

The CP violating 4 fermion operators are in- 
duced by neutral Higgs exchange at tree level [2, 
56] and, analogously to the 2-loop Barr-Zee con- 
tributions, they are sensitive to both the phases 
in the non-holomorphic Higgs couplings and the 
scalar-pseudoscalar mixing in the Higgs sector. 

Finally, contributions to the Weinberg 3 
gluon operator can be induced by 2-loop 
diagrams that are sensitive to the scalar- 
pseudoscalar mixing in the Higgs sector and 
in particular also to the phases in the non- 
holomorphic Higgs couplings [3, 57]. These 
contributions can be sizable, but they can be 
avoided to a large extent if the Higgs-top cou- 
plings are suppressed by non-holomorphic cor- 
rections. Keeping also in mind the large un- 
certainties in estimating the contribution of the 
3 gluon operator to the neutron EDM, we do 
not include it in our numerical analysis. Fol- 
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lowing this approach, we find that the most im- 
portant contributions are typically 2-loop Barr- 
Zee contributions to the mercury EDM and 
1-loop Higgsino-Wino-slepton contributions to 
the Thallium EDM in the regions of parameter 
space that we consider below. 

(iv) Constraints from direct Higgs searches at 
LEP, Tevatron and LHC are implemented using 
HiggsBounds [58] as well as the latest updates 
from Atlas and CMS [59]. Generically direct 
SM Higgs searches do not lead to strong con- 
straints, as the lightest Higgs boson is usually 
in the range 120— 140 GeV in the regions of pa- 
rameter space considered below. The two heav- 
ier Higgs bosons can be much heavier and, due 
to the moderate values of tan (3, SUSY Higgs 
searches are also not constraining possible large 
effects in B mixing. 

(v) Flavor observables: The main constraints 
come from AM<j and AM S as well as the branch- 
ing ratios of the decays B — > X s j and B s — > 

We also implement the constraints from 
B — > tv, B — > Dtv and K — > fiv. 

In Fig. 2 we show in two representative sce- 
narios the possible values of the NP phase in B s 
mixing in the M H + - tan j3 plane together with 
the above mentioned constraints. The plot on 
the left of Fig. 2 shows scenario I, where we 
chose a common sfermion mass of fh = 1 TeV 
and trilinear couplings A t = fh, A(, = A T = 0. 
We assume the absence of the higher dimcn- 



3 We remark that in a MFV framework, both the squark 
soft masses and the trilinear couplings can in general 
be expanded in powers of the flavor violating spuri- 
ons y q = Yq. Higher order terms in the expansion 
can in principle induce additional flavor changing ef- 
fects at the loop level. In the case of neutral meson 
mixing and B s — > /i + ^t~ however, the loop level ef- 
fects induced by the higher order tems will be sub- 
dominant compared to the dominant tree level contri- 
butions from the modified Higgs sector. Among the 
considered flavor observables, only the loop induced 
B — > X a -y decay can be noticeably affected by pos- 
sible higher order terms in the trilinear parameters. 
Considering such higher order terms would therefore 
add more flexibility in controlling the B —> X s ^f con- 
straint. As this constraint turns out to have only a 



sional operators (24) in the superpotential and 
set 

e Q = -0.8 , ei = 0, e 2 = 0.3 , 

e 3 = , e 4 = -0.5 - 0.8i , e t = 1 , 

a = lu = , 

fi = 200 GeV , m s = 1 TcV , 
M = 10 TeV . (42) 

Correspondingly, NP contribution to the B mix- 
ing amplitudes are generated through the Wil- 
son coefficient C4 and therefore effects are much 
larger in B s than in B^ mixing. We observe 
that the B s mixing phase can easily reach the 
2cr bound given in (18), 4^ p > —0.4, even for 
moderate values of tan /3 ~ 10 and very large 
Higgs masses of Mjj± > 1 TeV. In this region 
of parameter space, Higgs contributions to the 
B s — > fi + i-i~ decay and to B — > X s j are well 
under control. The strongest EDM constraints 
in this scenario come from the mercury EDM 
but due to the large Higgs masses they turn out 
to be easily fulfilled. As the higher dimensional 
operators in the superpotential are not present, 
vacuum stability bounds as well as electroweak 
precision constraints are always fulfilled. 

The plot on the right of Fig. 2 shows sce- 
nario II, where we allow both for higher di- 
mensional operators in the superpotential and 
Kahler potential, but consider new sources of 
CP violation only in the superpotential. We 
chose third generation squark soft masses of 
rriq L = mj B = mg H = 700 GeV, all remaining 
sfermion masses m — 1rriq L and trilinear cou- 
plings A t = 2 TeV, A b = A T = 0. In addition 
we set 

e = , ei=0, e 2 = -1.8 , 
e 3 = -1.6 , e 4 = l-6, e e = 1 , 
a = 1.5 — i , u = 1.5 , 
fi = 150 GeV , m s = 500 GeV , 



small impact on our analysis, the higher order terms 
would not change any of our conclusions. 
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FIG. 2. Possible values for the NP phase in the M H ± - tan (5 plane in the two example scenarios 
described in the text. The gray regions are excluded by AM S (solid line), AMd (long dashed line), BR(5 — > 
X s ^f) (dotted line), BR(B S — > fi + /i~) (short dashed line) and the mercury EDM (dash-dotted line). In the 
red hatched region of the right plot, the electroweak vacuum is not absolutely stable. 



M = 5 TeV . (43) 

In this setup the NP phases in B s and B^ mix- 
ing are induced by the Wilson coefficient C*2. 
Therefore they are of comparable size in both 
cases and typically well within the 2a ranges 
of (17) and (18). We observe that sizable val- 
ues for the B s mixing phase are possible even for 
small tan f3 ~ 4 but require a rather light Higgs 
spectrum, which agrees with our expectation 
that C2 can be important only for small Higgs 
masses. Because of the small values of tan/3, 
both constraints from B s — » H + H~ an d from 
the EDMs (that mainly come from the Thal- 
lium EDM) are well under control in the con- 
sidered scenario. Concerning the bound from 
BK(B — > X s <y) we remark that for the light 
Higgs masses, there are sizable charged Higgs 
loop contributions to the b — > amplitude 
that are further enhanced by higher order tan (3 
resummation factors. These contributions can 
be partly canceled either by tprns/M correc- 
tions to the couplings of the charged Higgs to 
the right-handed top quark, or by chargino-stop 



loops as long as stops are rather light, below 
1 TeV. In contrast to scenario I, a very impor- 
tant constraint is now coming from AM^. Also 
bounds from vacuum stability start to be im- 
portant. Contrary to the framework discussed 
in [14] however, vacuum stability bounds can be 
compatible with a large B s mixing phase with- 
out the need of additional physics that stabilizes 
the electroweak vacuum. 4 

In Fig. 3 we present the results of a parameter 
scan of the model. The left plot shows the cor- 
relation between the semi-leptonic asymmetries 
ctg L and ag L , while the right plot shows the cor- 
relation between the time dependent CP asym- 
metries S^Ks an d S^p^. Blue (dark gray) points 
correspond to a scenario where new sources 



4 We note that in the framework considered in [14] a 
large negative /i term is required to generate sizable 
flavor changing neutral Higgs couplings at the loop 
level. In the scenario considered here however, small 
values for fi are possible that considerably soften the 
vacuum stability bounds. 
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FIG. 3. Correlation between the semi-leptonic asymmetries ag L and ag L (left) and the time-dependent 
CP asymmetries S^k s and S^,^ (right). Blue (dark gray) points correspond to a scenario with new sources 
of CP violation in the Yukawa couplings, while green (light gray) points correspond to new CP phases in 
the Higgs sector. The red thin contours show the allowed regions using the bounds on (j)^ p and <^ p from 
Fig. 1. In the left plot, the magenta band (orange ellipse) shows the measurement of Ag h (ag L and a| L ) by 
DO. 



of CP violation are arising entirely from the 
modified Yukawa couplings and the higher di- 
mensional operators in the superpotential (24) 
are switched off completely. We consider mass 
scales as in scenario I above and allow tan/3 < 
15 as well as epsilon parameters < 1 with £o 
and eg real and arbitrary phases for the remain- 
ing £j. As expected, in such a setup the NP 
effects in B s mixing are much larger than in Bd 
mixing and the LHCb bounds on 4> s exclude a 
sizable NP phase in Bd mixing. 

Green (light gray) points correspond to a sce- 
nario where also the higher dimensional opera- 
tors in the superpotential (24) are considered. 
We fix mass scales as in scenario II above and 
allow the parameters \a\, \u>\ < 2 with arbitrary 
phases and real epsilon parameters |£j| < 2 as 
well as tan/3 < 15. In this setup, the NP con- 
tributions to CP violation in Bd and B s mix- 
ing are comparable. Consequently, a sizable 
NP phase in Bd mixing can be compatible with 
the LHCb constraints on the NP phase in B s 



mixing. Even though new CP phases come 
entirely from the superpotential operators, we 
stress that the presence of the Kahler poten- 
tial operators is crucial. They allow for low 
tan P values and therefore the constraint from 
B s — > fi + jU~ can be avoided. 

We conclude that in order to generate siz- 
able corrections to CP violation in Bd mixing 
that are in agreement with the LHCb data on 
CP violation in B s mixing, both the higher di- 
mensional operators in the superpotential (that 
modify the Higgs spectrum) as well as the 
higher dimensional operators in the Kahler po- 
tential (that lead to non-holomorphic Higgs- 
fcrmion couplings at tree level) are required. 

CONCLUSIONS 

Recent results from LHCb on the time- 
dependent CP asymmetries in B s — > tjxp and 
B s — > ipfo significantly restrict the allowed val- 
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ues for the B s mixing phase. Combining these 
results with measurements of CDF and DO of 
the time-dependent CP asymmetry in B s — > 
as well as the measurements of the time- 
dependent CP asymmetry in Bd — > ipKs at 
the B factories, we find the following 2er ranges 
for possible NP phases in B s and Bd mixing: 
-0.28 < <^ p < and -0.40 < < 0.20. 

The preference for a negative NP phase in Bd 
mixing is driven by tensions in fits of the Uni- 
tarity Triangle, while the NP phase in B s mix- 
ing is perfectly consistent with 0. Under the 
assumptions that neither the absorptive part of 
the B s mixing amplitude nor the SM tree-level 
B s iptfi decay amplitude are significantly af- 
fected by NP, the anomalous like sign dimuon 
charge asymmetry observed by DO cannot be 
explained given the above bounds. 5 

In view of these results we studied the possi- 
ble impact of higher dimensional operators in 
the MSSM on B physics. We considered di- 
mension 5 operators both in the superpotential 
and in the Kahler potential assuming that they 
are generated at a scale of M ~ 5 — 10 TeV. 
The 1/M suppressed operators in the super- 
potential have important impact on the Higgs 
spectrum. They can significantly enhance the 
tree level mass of the lightest Higgs boson and 
lead to a mass splitting between the two heavy 
neutral Higgs bosons. With complex coeffi- 
cients they also lead to scalar-pseudoscalar mix- 
ing. The 1/M operators in the Kahler po- 
tential can induce non-holomorphic Higgs cou- 
plings and consequently flavor changing neutral 
Higgs couplings at the tree level. Assuming 
that the Kahler potential operators follow the 
minimal flavor violation ansatz, we find a fla- 
vor phenomenology that resembles to a large 
extent the 2 Higgs doublet model with MFV 
discussed in [1, 2]. In particular large values for 



5 The possibility to reconcile the LHCb data on the 
time dependent CP asymmetries in B 3 — > rp<j> and 
B B — > ipfo and the DO measurement of the like sign 
dimuon charge asymmetry with NP effects in Ti2 has 
been analyzed very recently in [60]. 



tan /3 are not required to generate sizable correc- 
tions to B meson mixing from neutral Higgs ex- 
change and therefore the strong constraint from 
BR(£> S — » fi + n~) can be significantly relaxed. 

We find that both superpotential and Kahler 
potential operators are required to generate 
non-standard effects in the Bd mixing phase 
that are in agreement with the current bounds 
on the B s mixing phase from LHCb. The cor- 
responding region of parameter space is charac- 
terized by rather small values of tan j3 ~ 5 as 
well as low masses for the heavy Higgs bosons 
M H2 ,M Hs ~ 300 - 400 GeV. We stress that 
in the class of models discussed in this work, a 
non-standard Bd mixing phase does imply also 
non-standard effects in the B s mixing phase at 
a level that can be tested in the near future by 
LHCb. 
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APPENDIX 

Here we briefly present UV completions that 
lead to the higher dimensional operators in (24) 
and (25). The dimension 5 superpotential op- 
erator in (24) can be generated by integrating 
out a heavy singlet S with the following super- 
potential interactions [15] 

CD J d 2 6 + V^SH u H d ^J . (44) 

Additional gauge interactions that are broken 
at a high scale M can be effectively described 
by dimension 6 operators [15, 40]. 

One way to generate the Kahler potential op- 
erators in (25) is to introduce two heavy SU(2) 
doublets H u , Hd with hypercharge +1 and — 1 
that couple to fermions and mix with the MSSM 



13 



Higgs doublets H Ul H d [16]. 6 Neglecting for 
simplicity gauge interactions, their Kahler po- 
tential and superpotential read 

cd Jd^e (H\H d + hIhJ) 

+ J d 4 9 (H\H d + HlH u + h.c.) 
+ J d 2 e(MH u H d + \ t H d LE 

+ \ U H U QU + X d H d QD + h.c.) . (45) 

Integrating out the heavy Higgs doublets gener- 
ates the supersymmetric term in (25) and SUSY 
breaking can be incorporated with the auxiliary 
spurion Z . 

In addition to (25), integrating out the heavy 
Higgses from (45) also generates 1/M sup- 
pressed terms in the superpotential 

+ X u X e {QU)(LE)+h.c?j , (46) 

that in turn generate dimension 5 fermion - 
sfermion interactions. Such interactions are 
however not relevant for the topics discussed in 
the present work. 
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